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Abstract 

I discuss the theoretical uncertainties in the indirect Higgs mass determina- 
tion. I present the probabihty density function for the Higgs mass obtained 
combining together the information from precision measurements with the 
results from the direct search experiments carried out at LEP. The proba- 
bility that the Higgs weights less than 116 GeV comes out to be around 35% 
while the 95% upper limit is located around 210-230 GeV. 



*Talk presented at "50 Years of Electroweak Physics: A symposium in honor of Pro- 
fessor Alberto Sirlin's 70th Birthday", New York University, New York, October 27-28, 
2000 



1 Introduction 



The last months of the year 2000 have seen a great excitement in the physics 
community because of a possible evidence at LEP of a Higgs boson with a 
mass Mh 115 GeV. Unfortunately the shutdown of LEP will not allow 
additional information on the Higgs to be collected in the near future. Thus, 
it seems a good moment to try to review what we (do not) know about the 
Higgs. 

On one side we have the impressive amount of data collected at LEP, 
SLC, and the Tevatron that allow to probe the quantum structure of the 
Standard Model (SM), thereby providing indirect information about the 
Higgs mass. On the other we have the result of the direct searches performed 
at LEP with the excess of events reported in the combined analysis of the 
four collaborations. The virtual Higgs effects are usually analyzed through 
a fit to the various precision observables that allows a 95% Confidence 
Level (C.L.) upper bound to be derived. The outcome of the searches used 
to be reported as a combined 95% C.L. lower bound, while today is more 
appropriate described by the likelihood of the experiments. Given these two 
pieces of information it is often said that one of the greatest achievement of 
LEP has been to have pin down the Higgs mass between 113 (from the 95% 
C.L. lower bound of the direct searches) and 170 GeV (from the 95% C.L. 
upper bound of the global fit to clcctrowcak data). What I would like to 
discuss in this talk is how confident we are in the fact that 113 < < 170 
GeV addressing two aspects of this problem: i) how solid is the 170 GeV 95% 
C.L. upper bound? ii) What is the right way to combine the information 
from the precision measurements with that from the searches to answer the 
simple question: "what is the probability that the Higgs mass is between, for 
example, 113 and 170 GeV?" 

2 Theoretical uncertainties in the indirect Higgs 
mass determination 

As well known, the vector boson self-energies are sensitive to virtual Higgs 
effects and therefore the value of the Higgs mass affects the precision mea- 
surements like, for example, the ones made at peak of the Z. However, the 
Higgs behavior in the relevant electroweak corrections is very mild, just log- 
arithmic. To appreciate how much this logarithmic behavior makes hard to 
constrain the Higgs let me consider a simple example. Let us take the effec- 
tive electroweak mixing angle, svcP'O^^J^^ = Sg//' ^^^^ most important 
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quantity in the determination of M^, and write it as 



„2 



(ci + 6ci) + (C2 + 5c2)log y- y = (M^/lOOGeV). 



(1) 



In Eq. (|T|) I identify the l.h.s. with the experimental result s 



eff 



eff)0 



a{s^jj) while in the r.h.s. 5ci represent the theoretical uncertainty in the 
corresponding coefficients connected to the fact that we have computed Cj 
in perturbation theory through certain order in the perturbative series and 
therefore we do not know their exact values because of higher order contri- 
butions. From Eq. one obtains 



y = 2/0 exp 



C2 C2. 



Ath = 5ci + 5c2 log y 



where yo is the value corresponding to Sci = 5c2 
of Af/j in extracting Mh I put cr = and take 



(T = 0. To see the effect 
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5.5 X 10^ 



Ath ~ ±1-4 X 10" 



(3) 



where ~ 0.23, (? = 1 — s^. In Eq. (|3|) I estimate C2 through the Higgs 
leading behavior of the correction Ar relevant for Sg^j [|l| while for Ath, I 
take the value estimated in the 1995 CERN report on 'Precision calculation 
for the Z resonance' H] that was supposed to represent the uncertainty due to 
next-to-leading two- loop electroweak effects. Inserting the values of Eq. (^) 
into Eq. @ yields y ~ 1.29 yo- We see that a theoretical uncertainty coming 
from two-loop unknown contributions (that are supposed to be not even the 
dominant part) makes an error in the indirect determination of the Higgs 
mass of 29 %! 

Eq.(||) tells us an important thing, namely that the error on y depends on 
its central value (yo)- This implies that not always shrinking the uncertainty 
in Iny reduces the uncertainty on y. This is true only if the central value 
does not change, but this is not always the case (as we will see later). This 
consideration can be put in a more formal way noticing that if the logarithm 
of a quantity, A = \ny, is normally distributed, then the quantity itself is 
distributed according to a lognormal (see, e.g., Q for the properties of this 
distribution) whose standard deviation, given by 



a{y) = (exp [2E[A] +2a'^{A)\ - exp [2E[A] + a'^{A) 



is a combination of the expected value and standard deviation of its loga- 
rithm and therefore compensating effects can happen. 
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The above example clearly tells us that to extract accurate indirect in- 
formation on the Higgs one needs both very precise experiments and a very 
good control of the theory side. This brings in the issue of what error we can 
associate to our theoretical predictions. They are affected by uncertainties 
coming from two different sources: one that is called parametric and it is con- 
nected to the error in the experimental inputs used in our predictions. The 
second one is called intrinsic and it is related to the fact that our knowledge 
of the perturbative series is always limited, usually to the first few terms. 
Concerning parametric uncertainties, a(0), and Mz are very well mea- 
sured, while the top mass, M^, and the strong coupling constant, a^, are not 
so precisely known. However, the scale of the weak interactions is given by 
the mass of the intermediate vector bosons, so what actually matters in our 
predictions is not a(0) but a{Mz)- The latter contains the hadronic contri- 
bution to the photon vacuum polarization, {Aa)h, that cannot be evaluated 
in perturbation theory. Fortunately, one can use a dispersion relation to 
relate it to the experimental data on the cross section for e^e~ annihilation 
into hadrons. In the recent years this subject has received a lot of attention 
with the appearance of several new evaluations that have followed two main 
streams. On one path there are the most phenomenological (ph.) analyses 
that rely on the use of all the available experimental data on the hadron 
production in e"*" e~ annihilation and on perturbative QCD (pQCD) for the 
high energy tail {E > 40 GeV) of the dispersion integral 0, || . On the other 
the so called "theory driven" (t.d.) analyses that advocate the use of pQCD 
down to energy scale of the order of tau mass, supplemented by the use of 
the experimental data in regions, like, for example, the threshold for the 
charmed mesons, where pQCD is not applicable ^, §]• I am not going to 
discuss the differences in the various analyses (see Fred Jegerlehner's talk 
[^) but I would like to point out few facts: i) all results are compatible with 
each other. The choice of one value instead of another is just a matter of 
taste (or friendship), ii) The t.d. results have a smaller error with respect 
to the ph. ones but a lower central value, iii) I can take two perfectly com- 
patible numbers, one from the ph. analyses like {Aa)h = 0.02804 it 0.00065 
1] and one from the t.d. ones like (Aa),, = 0.02761 ± 0.00022 @ and get 
a difference in the 95% C.L. upper bound for the Higgs mass 0(50 GeV). 
Just to mention, it is the t.d. value, that has a smaller error, that gives the 
higher upper bound. 

I would like to emphasize that this uncertainty has nothing to do with 
the blue band in the famous Ax^ vs. LEPEWWG plot. There, the 
blue band represents, for a chosen value of {Aa)h, the intrinsic uncertainty. 
There is no way to rigorously define the intrinsic uncertainty. The best that 
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can be done is to try to estimate it by comparing the output of the two codes 
TOPAZO and ZFITTER that now include up to two-loop next-to- 
leading terms, when they are run enforcing the several built in options for 
resumming known effects. These options are supposed to mimic the size of 
unknown higher order terms and the numerical spread in the outputs of the 
two codes can be taken just as an order of magnitude of the unknown higher 
order contributions. 



3 Higgs mass inference from precision measure- 
ments 

I am going to discuss now what can we learn about the Higgs from precision 
measurements. In the spirit of the question raised in the Introduction, I 
am not going to perform the standard analysis, although clearly what I 
present will be related to it, but following a Bayesian approach I am going 
to construct /{mn \ ind.), the p.d.f. of the Higgs mass conditioned by this 
indirect information under the assumption of the validity of the S.M. I would 
do it employing the three observables, -Sejj, M^^, and F^. These quantities 
are the most sensitive to the Higgs mass and also very accurate measured. 
The most convenient way to approach the problem is to make use of the 



simple parameterization proposed in Ref.|12| and updated in Ref.|13] where 



Sgjj, My^r and Ti are written as functions of Mjj, Mt, Og and the hadronic 
contribution to the running of the electromagnetic coupling: 

slff = {sljf)o+ClAi + C2A2- C3A3 + C4A4, (4) 

Mw = M^- diAi - d^Aj - d2A2 + ds^s - d^A^ , (5) 
= - gi Ai - 95 - 92 A2 +gsA3- A^ . (6) 

In the above equations Ai = ln(MH/100GeV), A2 = [(Aa)/,/0.0280 - 1], 
^3 = [(Mj/175GeV)2-l] and ^4 = [(a,(Mz)/0.118 - 1]; (s2^^)o, 
and are (to excellent approximation) the theoretical results obtained at 
the reference point {Aa)h = 0.0280, Mt = 175 GeV, and a^(Mz) = 0.118 
while the values of the coefficients Cj, di and g^ are reported in Tables 3-5 of 



Ref.[13| for three different renormalization scheme. Formulae (|^-|6[) are very 
accurate for 75 <^ Mjj <^ 350 GeV with the other parameters in the ranges 
170 ^ Mt <, 181 GeV, 0.0273 <, (Aa)^ ^ 0.0287, 0.113 <, a,{Mz) <, 0.123. 
In this case they reproduce the exact results of the calculations of Refs.([ll, 



15L 13(1) with maximal errors of Sslfj ~ 1 x 10-^ 5M„ <:, 1 MeV and 6Tp <:, 3 



KeV, which are all very much below the experimental accuracy. Outside the 
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above range, the deviations increase but remain very small for larger Higgs 
mass, reaching about 3 x 10"^, 3 MeV, and 4 KeV at Mh = 600 GeV for 
slj^j, Mw, Te, respectively. 

Formula (§) can be seen as providing an indirect measurement of Ai = 
[s^jj: — (Se//)o ~ C2A2 + C3^3 — CiA^l/ci, while (|5|) and (^) of the quantities 
Y = M^--M^- d2^2 + ^3^3 - diAi and Z = r° - - 52^2 + 53^3 - 54^4 , 
respectively, all three variables being described by Gaussian p.d.f.'s. The 
Ai, Y and Z determination are clearly correlated, therefore one has to built 
a covariance matrix. This can be easily done because formulae (0-|6|) are 
linear in the common terms X = {Mf, as{Mz), {Aa)h}- The likelihood of 
these indirect measurements = {Ai, Y, Z} is then a three dimensional 
correlated normal with covariance matrix 

or 

/(^l ln(mH/100)) oc e-^'/2 (g) 

where = A^V^^A with A^ = {ai - ln(m///100), y - di ln(mH/100) - 
d5ln^(m///100), z - 51 ln(m///100) - g^ln^imn /WO)}. Using Bayes' theo- 
rem the likelihood can be turned into a p.d.f. through the choice of a 
prior. The natural choice is a uniform prior in ln(m_H'), since it is well under- 
stood that radiative corrections measure this quantity (for this reason the 
likelihood (|8|) has been expressed in terms of ln(?n,//)). The uniform prior in 
ln(m//) implies that /(ln(m/f/100) | ind.) is just the normalized likelihood 
(P). Thus, using standard probability calculus I can express the result as a 
p.d.f. of Mh: 

/{iJiH I ind.) = „^ -(y2/2)A • 

Jo ^ dniH 

4 Including the constraint from the direct search 

Given /{uih \ ind.) a natural question to ask is how this p.d.f. should be 
modified in order to take into account the knowledge that the Higgs boson 
has not been observed (that there is some indication for the Higgs) at LEP0. 
To answer this question let me discuss first an ideal case. I consider a search 
for Higgs production in association with a particle of negligible width in 
an experimental situation of "infinite" luminosity, perfect efficiency and no 

related discussion can be found in Ref. 1191. 
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background whose outcome was no candidate. In this situation we are sure 
that all mass values below a sharp kinematical limit Mk are excluded. This 
implies that: a) the p.d.f. for Mh must vanish below Mk] b) above Mk 
the relative probabilities cannot change, because there is no sensitivity in 
this region, and then the experimental results cannot give information over 
there. For example, if M^- is 110 GeV, then /(200 GeV)//(120 GeV) must 
remain constant before and after the new piece of information is included. 
In this ideal case we have then 

r niH <Mk 

fimH\dtr.kind.) = l nin > Mk , ^^^^ 

L JMif /("iff I ina.) CLmH 

where the integral at denominator is just a normalization coefficient. 

More formally, this result can be obtained making explicit use of the 
Bayes' theorem. Applied to our problem, the theorem can be expressed as 
follows (apart from a normalization constant): 

/{niH I dir. & ind.) oc f{dir. \ mn) ■ f{mH \ ind.) , (11) 

where f[dir \ m//) is the so called likelihood. In the idealized example we are 
considering now, f{dir \ rrifj) can be expressed in terms of the probability of 
observing zero candidates in an experiment sensitive up to a Mk mass for 
a given value rriH, or 

fidzr. I mn) = /("zero cand." | m^) = | ° t^/ (^2) 

11 niH > Mk ■ 

In fact, we would expect an "infinite" number of events if Mh were below 
the kinematical limit. Therefore the probability of observing nothing should 
be zero. Instead, for Mh above Mk, the condition of vanishing production 
cross section and no background can only yield no candidates. 

Consider now a real life situation. In this case the transition between 
Higgs mass values which are impossible to those which are possible is not 
so sharp. In fact because of physical reasons (such as threshold effects and 
background) and experimental reasons (such as luminosity and efficiency) 
we cannot be really sure about excluding values close to the kinematical 
limit, nevertheless the ones very far from Mk are ruled out. Furthermore, 
the kinematical limit is in general not sharp; at LEP, for example, the large 
total width of the Z° plays an important role. Thus, in a real life situation 
we expect the ideal step function likelihood of Eq. (|12|) to be replaced by 
a smooth curve which goes to zero for low masses. Concerning, instead, 
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the region of no experimental sensitivity, Mjj ^ Mx^ff, the hkehhood is 
expected to go to a value independent on the Higgs mass that however is 
different from that of the ideal case, i.e. 1, because of the presence of the 
background. 

In order to combine the various pieces of information easily it is con- 
venient to replace the likelihood by a function that goes to 1 where the 
experimental sensitivity is lost Because constant factors do not 

play any role in the Bayes' theorem this can be achieved by dividing the 
likelihood by its value calculated for very large Higgs mass values where no 
signal is expected, i.e. the case of pure background. This likelihood ratio, 
TZ, can be seen as the counterpart, in the case of a real experiment, of the 
step function of Eq. (12). Therefore, the Higgs mass p.d.f. that takes into 
account both direct search and precision measurement results can be written 
as 

ff I „ . ,x n{mH) f{mH\ind.) 

f{mH dir. kind.) = — . 13 

Jo Tl{mH) J [rriH \ md.) dmn 

In Eq. ( p!3|) TZ, namely the information from the direct searches, acts as a 
shape distortion function of /{rriH \ ind.). As long as TZ{mH) is 1, the shape 
(and therefore the relative probabilities in that region) remains unchanged, 
while TZ{mH) — > indicates regions where the p.d.f. should vanish. 

One should notice that TZ{mH) can also assume values larger than 1 for 
Higgs mass values below the kinematical limit. This situation corresponds to 
a number of observed candidate events larger than the expected background. 
In this case the role played by TZ{mH) is to stretch f{mH \ ind.) below the 
effective kinematical limit and this might even prompt a claim for a discovery 
if TZ becomes sufficiently large for the probability of Mh in that region to 
get very close to 1. 



5 Results 



The experimental inputs I use to construct f{mH \ ind.) are [1£]: s^^j 



0.23146 ± 0.00017, = 80.419 ± 0.038 GeV, = 83.99 ± 0.10 MeV, 
Mt = 174.3 ± 5.9 GeV, as{Mz) = 0.119 ± 0.003. Concerning (Aa)/,, as I 
said we have many possible choices. I am going to present my results using 
two different values of {^a)h-, one as representative of the ph. analyses and 
the other for the t.d. evaluations. For the ph. case I take the standard value 
(Aa)f = 0.02804 ± 0.00065 while for the t.d. analyses I choose the one 
that has the smallest uncertainty, i.e. (Aa)^^ = 0.02770 ±0.00016 0. The 
intrinsic uncertainty is taken into account by averaging different inferences 
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Table 1: Summary of the direct plus indirect information, 
(see Ref.0). 

The values of the TZ function that enters in Eq. (p^) has been provided 
by the LEP Higgs Working Group [^0|: they take into account all Higgs 
searches by the four LEP collaborations. 

Table I summarizes the result of my analysis. The shape of the p.d.f. 
with and without the inclusion of the direct search information is presented 
in Fig. |l[ From this figure one can notice that, in the case of f{niH \ ind.), 
the use of a higher central value for (Aa)^ (i.e. (Aa)^"^) tends to concentrate 
more the probability towards smaller values of Mh- As a consequence the 
analysis based on {Aa)j^'^ gives results for the expected value, E[M^], the 
standard deviation, a{M„), of the Mh p.d.f. and 95% probability upper 
limit, Mfj, very close to those obtained using (Aa)^^ regardless the fact 
that the uncertainty on (Aa)^"^ is approximately 4 times larger than that 
of (Aa)^^. As expected, the inclusion of the direct search information 
drifts the p.d.f. towards higher values of Mh by cutting regions below 110 
GeV. However, the most striking thing in the figure is the spike around 
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100 150 200 Mh 




100 150 200 



Mh 



Figure 1: Probability distribution functions using only indirect information 
(dashed line) and employing also the experimental results from direct searches 
(solid one): a) {Aa)h = 0.02804(65); b) {Aa)h = 0.02770(16). 



M[{ « 115 GeV. Indeed, due to the excess of events in the combined data 
set of the LEP experiments, the TZ function around 115 GeV takes values 
~ 30 strongly enhancing the probability of that region. 

6 Conclusions 

The analysis that I have presented clearly shows that a heavy Higgs scenario 
is highly disfavored, given a 0(90%) probability that the Higgs weights less 
than 200 GeV. However, it should be said that the distribution has a quite 
long tail. Indeed we have that above 300 GeV there is still a residual C(l%) 
probability. The excess in the data recorded by the LEP collaborations 
is reflected in the analysis by the spike in the distribution around ^ 
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115 GeV. Maybe it can be of some interest to know that, according to my 
analysis, the probabihty that the Higgs mass is below 116 GeV is 37%, if we 
use (Aa)f and 38% in the case of {Aa)j^^ . 

Let me say clearly that all the results I have presented are derived under 
the assumptions of the validity of the SM and rely on the experimental inputs 
I have used. In particular I have Tiscd for s^jj the combined LEP+SLD value, 
although it is well known that the two most precise determinations of it are 
not in very good agreement. Just to see how much my conclusions rely on the 
combined s^jr^. value I have redone the analysis discarding the SLD result, 
i.e. using the LEP value for the effective sine, s^jf = 0.23184 ± 0.00023. 
The results of this extreme exercise are not very different: in the worst case, 
i.e. employing (Aq;)^^, P{Mh < 200 GeV) = 75% while Mf = 320 GeV. 
Finally there is still a 23% probability that the Higgs mass is below 116 
GeV. 



I am grateful to G. D'Agostini, S. Fanchiotti, F. Feruglio, P. Gambino, 
M. Passera, A. Sirlin and A. Vicini for the fruitful collaboration on different 
projects variously connected with the subject discussed here. I would like 
also to thank M. Porrati and the other organizers for the excellent organi- 
zation and the very pleasant atmosphere of the workshop. This work was 
supported in part by the European Union under the contract HPRN-CT- 
2000-00149. 
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